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Abstract 

In the theory of gelation it has been implicitly assumed that (I) a cyclic bond is a finite bond that returns to itself; 
(II) cyclic bonds distribute at random in network structures. In this paper these two assumptions are reexamined 
from a new point of view. The physical soundness of the assumptions is assessed through comparison with experi¬ 
mental observations. 
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In developing the theory of branching processes |U[2l[4Mj, we have introduced two basic assumptions: (I) a cyclic 
bond is a finite link and (II) cyclic bonds distribute randomly in network structures. These two assumptions have 
served as the basis of the theory of gelation, leading us to the understanding of various phenomena that occur in 
the gelling processes: upward shift of gel points, existence of critical dilution, and occurrence of permanent sol. In 
this paper, we reinvestigate the above basic assumptions in order to lay those on firm foundation. 

1 Theoretical Background 

1.1 On Assumption I 

In the theory of branching process, it has been implicitly assumed that a cyclic bond is a finite link, because it 
returns to the starting point. To date this statement has been accepted as a theorem. However, when analyzed 
in detail, it is seen that the statement is not necessarily self-evident. To make this issue clearer, an example of 
a branched structure involving one cyclic bond is illustrated in Fig. [l] for the multiple link system of / = 3 and 
J = 3: large filled circles (O) represent monomer units, open circles (o) functional units and the symbols (cQ) 
junction points. Suppose that a cyclization has occurred to create the bond 1 — 6. Taking a look at this figure, one 
may raise the objection: 

Why can the author affirm that a cyclic bond is finite? It is certain that the bond 1 — 6 emanating from the 
junction point 1 returns to itself through the route 1 — 6 — 2 — a—1, but at the same time it can extend to infinity, 
for instance, through the route 1 — 6 — 3, doesn’t it? 

That’s right! But the true meaning of the cyclic bond resides in a deeper level. To discuss this problem, it is 
necessary to keep in mind the following two facts: (A) a cyclic bond gives no effect on the growth of cluster size; it 
only wastes functional units; (B) once a cyclic bond is formed, the memory of the event of cyclization is completely 
lost; namely, anyone of 4 bonds (1 — a, a — 2, 2 — 6 and 6—1) can be a cyclic bond. 

Because of the fact (A), there must be at least one useless, finite bond that returns to itself. Let it be the 
bond 1 — 6 as defined above. Then disconnect the bond 1 — 6, as shown in Fig. [jJ(B). The situation now becomes 
apparent. The path from the junction point 1 to 3 is still alive through another route 1 — a — 2 — 6 — 3. Whether the 
bond 1 — 6 is closed or open, the route to the junction point 3 is already constructed through the (intermolecular) 
bond 1 — a. It is seen that the cyclic bond is a bond not emanating from the junction point 1, but returning to 
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itself (see Fig. QJ(C)). When we state that the bond 1-b is a finite bond, it is a cyclic bond, whereas when we state 
that the bond 1-b can lead to infinity, it is an intermolecular bond. Both the statements are true because of the 
fact (B). The same argument, of course, applies equally to the other bonds, 1 — a, a — 2 and 2 — 6. 

1.2 On Assumption II 

As far as we consider individual clusters such as shown in Fig. [T] it is obvious that the validity of Assumption 
II is confined to the cyclic structure l-a-2-b-l. Then it might appear that Assumption II can not be extended 
beyond the local structure, since individual clusters should contain different sizes, different numbers and different 
types of rings, so the cyclic structures change from clusters to clusters. Now the problem with which we are faced 
turns out to be much intricate, theoretical treatment appearing harder. However, in most theoretical treatments 
of branching process, what is required is the mean properties of the system; cases as to require detailed structural 
informations of individual clusters are rare. In such average pictures [6], every functional unit has an equal chance 
of cyclization or intermolecular reaction. Thus, taking the above-mentioned fact (B) into consideration, it may be 
concluded that there is a sufficient reason for us to accept Assumption II, the random distribution of cyclic bonds. 




Fig. 1: Representation of a Cyclic Bond in the Multiple Link System of / = 3 and J = 3. 

Our task is then to make a thorough assessment of the validity of Assumptions I and II through comparison 
with extensive experimental observations. 

2 Formulation of Gelling Processes 

Consider the branching process in the mixing system of the R—A f + R—B g model comprised of various monomers 
having two different type of functional units, {/iMa*} and where Ma ; and Mb,, are the numbers of the 
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A and B type monomers, respectively; /,; and gj are the corresponding functionalities having positive integers with 
the subscripts, i and j, signifying 1,2,3, •••. Chemical bonds can be formed only between A and B functional 
units. Let [r] be the total molar concentration of cyclics and V the system volume. Bearing in mind that each 
ring possesses only one cyclic bond [Bj, we introduce the notations: 


• fraction of cyclic bonds as against A functional units : 

• fraction of cyclic bonds as against B functional units : 

• fraction of A,; functional units : 

• fraction of Bj functional units : 


pr. A = [r]/ E s f* M a.) = [r }/c f 

pr b = [r]/ (vtv E, 9 3 Mb 0 ) = [T}/C g 
XK =/ i M Ai /V./ i M Ai 

A - J t 

Mb =9jMB j /^2 j 9jM Bj 


( 1 ) 


2.1 Gel Point 

Given Assumptions I and II, the number, a, of branches leading from a given R— Af unit, via R— B g units, to the 
next R— Af units is simply given by 


<* = E**i E k ( 9 \ ^Pb^-Pb) 3 *-'-^ ~ Pb(Z\X)} 

j k =0 ' ' 

x E XA, E 4 fl 7 1 W 1 - PA)^- 1 -^! - Pa(Z\X)} (2) 

i e=o \ / 

where Pa{Z\X) = pr a /pa and P B (Z\X) = pr b /p b represent the conventional conditional probabilities. The 
critical condition occurs at a = 1, so that 

(pac - Pr Ac ) (pbc - Pr Bc ) = -rrr\ — ttt, —r —tt ( 3 ) 

((/w) - l)((5w) - 1) 

Let k = 9 jM Bj / J2i fiMAi be the molar ratio. Then eq. © may be recast in the form: 

PAC ^ «/w) - lK(ffw) - 1) + W (4) 

— PAcq + Pr Ac 


which is of the form: p = p(inter) + p(ring), as expected. Note that the p(ring) term in eq. 0 is the quantity at 
Pc so eq. 0 indicates that the gelation occurs at the point where the ratio of the intermolecular bonds to the 
total number of possible bonds attains the classic gel point with no rings. 

To solve eq. 0 we must express [T c ] as a function of pac- It is well-known that the concentration of cyclic 
species is held constant above the critical initial-monomer-concentration, Co > Cg [6]. For this reason we may 
regard the following liming solution of Co —> 00 as the general solution of [T]: 


[r] = E ip - 


[«/w)-i)((gw)-iK/4 

2N Av x 


(for Co > C*) 


(5) 


Unfortunately, as emphasized previously (6j, eq. (0 breaks down for pa > PAc 0 - So we cannot link directly the 
limiting solution 0 with eq. 0. To resolve this problem, we expand eq. © about pa = PAc 0 to obtain 


[T(pa)] * E 

X=1 


Vx 


2N Av x 


«/w) - l)((gw) - 1) ^ <p t 


^-( PA - PAco ) 
, ^Av 
x=l 


( 6 ) 


In eq. 0, we have made use of the classic relation: ((/ w ) — l)((g w ) — P)P 2 Aco /r = 1. Let s = ((/ w ) — l)((g w ) — l)//c. 
Then replace P a in eq. © with PAc and substitute the resulting equation into eq. 0 to yield 


PAc 



\ - 


V~s 


Nj 


Av 


EL Vx (1 ^ h) 7/ 


yfs 


(7) 
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where 7 / = 1/C/ is the reciprocal of the initial A functional unit concentration. If we define the total functional 
unit concentration Cf g = N ^ v ( Y/i fiM a 4 + Yhj 9 jMBj ) and the total monomer concentration C = N ^ v ( Y2i Ma, + 
Y2j Mbj), the following relationship holds between 7 /, 7 / s = 1 /C/ g and 7 = 1/C: 


7/ = (1 + k)7/ s = 


ifn)K+{gn) 
(fn){g n ) 7 


(l/mol) 


( 8 ) 


where the subscript n denotes the number average quantity. 

The condition of the critical dilution can be obtained by imposing the boundary condition, pac 0 < PAc < 1, on 
eq. 0 : 


0 < 7/ < 7/c 


_ 1-l/Vi 

Tib Xb + Vs + l/2x) 


(9) 


It is important to notice that eq. ([9]) is meaningful only if k > 1, since otherwise the boundary condition, pa c = 1, 
cannot be fulfilled. 


2.2 Post-gelation 

Let Qa be the probability that a chosen branch emanating from the R—A/ monomer unit is finite and Q b be the 
corresponding probability for the R—B g monomer unit. Given Assumptions I and II, one can readily describe a 
set of the recurrence relations: 

Qa =1-pa+Pa ^.Xb, \Pb(Z\X) + (1 - P b (Z|A))Q® 3 ' -1 ] 

J (10) 

Qb =1-Pb+Pb XA S [Pa(Z\X) + (1 - Pa^IA))^/*- 1 ] 

The solution other than Qa = 1 or Qb — 1 is 

P aPb(1 - Pa(Z\X))(1 - P B (Z\X)) xa, £ QaJ ^ . Xb.., £ Q^ = 1 (11) 

With P A (Z\X) = pr k /pa and P B (Z\X) = Pr b /pb, and hence pr a /pa = Pr b /pb in mind, substituting Q A = 
Qb = 1 into eq. CD, we recover the foregoing critical condition: 


PAc \l ((/ w ) - iK(sw) - 1 ) +PRac 

with k = gjM Bj / fi^Ai = Pa/pb as defined above. 

The sol and gel fractions of this system have the forms: 

W so i = . w Ai Qa + yb • WB i 

1 J 

W gel = 1 - W sol 


0 ) 


( 12 ) 


where WAi and w Bj represent the weight fraction of the ith A-type monomer unit and that of the jth B-type 
monomer unit, respectively: 


__ m Ai M Ai _ 

Li TO A i -Ma 4 + M Bj 

= __ 

’ 3 V, III A, A / A, + "'IS, A/b. 


(13) 


and satisfy w A t + Y2j w Bj = 1 - 


3 Examination of Results 

Eqs. (4), (7), (10) and (11) are the direct consequence of Assumptions I and II. In this paper we focus our attention 
to the post-gelation problem to reexamine the physical soundness of our assumptions through comparison with 
experiments. 
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For the assessment of the above theoretical reasonings, we take up the Ilavsky-Dusek observations mm for the 
gelling process of 4,4 / -diphenylmethane diisocyanate (DPMDI) and LHT240 at 55 ~ 60 °C in non-solvent state (see 
Fig. [2J|. 

Let DPMDI be the II—Aj unit (/ = 2) and LHT240 the R—B g unit, respectively, and we have re = [OH]/[NCO]. 
LHT240 is prepared from 1,2,6-hexanetriol by means of the anionic polymerization and is considered to be a 
mixture of diols (gi = 2) and triols (<? 2 = 3). Let x\ = M Bl / JT M B . be the mole fraction of the diols and 
x 2 = Mb 2 / Mbj that of the triols to the total amount of alcohols, so that x\ + X 2 = 1. According to the 
literature [1|, the mean functionality is (g n ) = 9j x j = 2.89. Thus x\ = 0.11 and x 2 = 0.89. These give 
(< 7 w) = ( g 2 )/{gn ) = 2.92. Clearly {xj} can be linked with the aforementioned quantities {xb,} by the equation: 

XB;, = 9jXj/^2.gjXj (14) 



h 2 c— o 

HC-O 

CH 2 

1 | J 3 

h 2 c— o 



CH 3 



H 

CH 3 



H 

C- 

CH 3 


O 

O 


O 



DPMDI 


-£=> 



arm 1 

arm 2 


arm 3 


Polyurethane Network 


Fig. 2: Formation of polyurethane network by the reaction of DPMDI and LHT240. 



H s O H ch 3 

Fig. 3: Backbone structure of the polymer from DPMDI and LHT240. 


For the present case, XBi = 2 • 0.11/2.89 = 0.076 and xb 2 = 3 • 0.89/2.89 = 0.924. Solving eqs. (fT 0 l) - (fTTl) under the 
condition mentioned above (/ = 2, g\ = 2, c / 2 = 3), one has 


together, by eq. ©, with 


Qa =1 


re {re - (I + XB 2 HPA -Pr a ) 2 } 


Q B = 


Xb 2 (pa — Pr a ) 3 
w - (pa-Pr a ) 2 
X'b 2 (pa — Pr a ) 2 


(15) 


F(pa)] I ^ 
PR A = — } ~ = \ 2^ 


Px 


2Na v x 


DU I 

+ \J ((/w) - l)((<?w) - l)/« 22 ~ PAc 0) f 1f 


^ ai—1 X— 1 

Using the above quantities, the gel fraction can be calculated by the equation 


(16) 


W ge i = 1 - ( W A Ql + W Bl Q£ + Wb 2 Qb) 


(17) 
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and 


w a = 


w Bi = 


\9n/'m A 


(gn)m A + fn{m Bn ) 
fnXjmBj 


(18) 


(. 9n)m A + / k ( to b „) 

For the diols (tobi) and the triols (m B 2 ) in question, the observed polydispersity index (m BtW ) /(m B , n ) = 1.03 [2j 
gives two possible solutions: (to B i,to B 2 ) = {(359, 751), (1057,665)}. In this paper we examine the solution, 
(m A , wbi, m B 2 ) = (168,1057,665). 

Before performing our simulation, we must evaluate the relative cyclization frequency tp x defined by 

nd/2v x 

^ x =(d/2n d ^) ti-'e-'dt (ig) 


Fz = i r l)/ l s = Op^eX- 1) 


For this purpose, it is necessary to acquire the information about the effective bond number £ e within the repeating 
unit. Set the standard bond length l s = 1.36 A, and the virtual bonds, I 3 = I 4 = 5.69 A (see Fig. Then we can 
calculate the required number for the (i,j) arm pairs as 


= —^{1.36 2 + 1.36 2 + 5.69 2 + 5.69 2 + 1.36 2 + 1.36 2 + (3.3)(1.41 2 + 1.53 2 + 1.41 2 ) 

+ 1.41 2 + 1 x 1.53 2 + 1.41 2 + (3.3)(1.41 2 + 1.53 2 + 1.41 2 )} (20) 

where t = 1 for £^ 2 , t = 4 for £ 2,3 and t = 5 for £ 1 ^ 3 . And we have 

fe = i(6.2+ 6,3 +6.3) = 67.92 (21) 

so that v x = (r 2 )/! 2 ~ Cf(6Sx — 1). We set CV = 4.5, the same value as that for the HMDI-LHT240 system [B . 
Now we can evaluate the relative cyclization frequency kp x under the condition of d = 3. The results are summarized 
in Table [Q Making use of these data, we can plot the weight fraction W ge i of gel as a function of the molar ratio 
k = [OH]/[NCO]; the result is illustrated in Fig. [^together with the experimental points by Ilavsky and Dusek [I]; 
the solid line represents the theoretical line based on eqs. GSD-dD, and the dotted line the prediction by the 
ideal tree theory without rings (pn = 0). Agreement between the theory and the experiment is vey excellent for 
k > 1, whereas marked disagreement is observed for the n < 1 zone. The latter phenomenon may be ascribed to 
side reactions, under the condition rich in —NCO moiety, caused by the recombination between urethane bonds 
—NH(CO)0— and unreacted —NCO functional units to give rise to, for instance, allophanate structures [Tj. The 
results are in good accord with the theoretical calculation based on the cascade theory by Ilavsky and Dusek [l]. 

It may be concluded that the present results confirm the physical soundness of Assumptions I and II. 


Table 1: Physicochemical parameters for DPMDI—LHT240. 


parameters 

unit 

values 

molecular mass 

m 

DPMDI : m A = 250 LHT240 : (m B ,„) = 708 

specific gravity (60 °C) 

p 

DPMDI : p A - 1.0 LHT240 : p B ~ 1.0 

/ 


2 

<0n> 


2.89 

<Sw> 


2.92 

c F 


4.5 

Ze 


68 

Is 

(A) 

1.36 

cyclization frequency 



E“=i 

(■ mol/l ) 

0.108 

00 1 

2^x=l Vx 2 N Av x 

(mo/ /l) 

0.028 
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Fig. 4: Gel fraction as against k = [OH]/[NCO]. Solid line ( —): theoretical line by eqs. (1151) — (1191) : dotted line (•••) : 
prediction by the ideal tree theory with no rings (pr = 0); (o): experimental points by Ilavsky and Dusek [lj. 
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